Texture and shape of two-dimensional domains of nematic liquid crystal 
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We present a generalized approach to compute the shape and internal structure of two-dimensional 
nematic domains. By using conformal mappings, we are able to compute the director field for a given 
domain shape that we choose from a rich class, which includes drops with large and small aspect 
ratios, and sharp domain tips as well as smooth ones. Results are assembled in a phase diagram that 
for given domain size, surface tension, anchoring strength, and elastic constant shows the transitions 
from a homogeneous to a bipolar director field, from circular to elongated droplets, and from sharp 
to smooth domain tips. We find a previously unaccounted regime, where the drop is nearly circular, 
the director field bipolar and the tip rounded. We also find that bicircular director fields, with foci 
that lie outside the domain, provide a remarkably accurate description of the optimal director field 
for a large range of values of the various shape parameters. 

PACS numbers: 68.55.-a, 61.30.-v, 68.60.-p 



I. INTRODUCTION 

Nematic liquid crystals are uniaxial fluids that re- 
spond elastically to an externally imposed deformation 
of the director, i.e., the locally preferred orientation of 
the molecules. The interfacial properties of nematic liq- 
uid crystals are also anisotropic. This means that the 
surface free energy depends on the attack angle of the 
director relative to the interface [THE]. Under conditions 
of spatial confinement, i.e., in a cavity surface and bulk 
elastic properties compete so as to minimize the overall 
free energy [9HTT]. 

For instance, homeotropic anchoring of the liquid crys- 
tal to the interface in a spherical drop or cavity may lead 
to a radial director field, whereas planar anchoring pro- 
duces a bipolar field EJJ . The former is characterized by 
a hedgehog point defect in the center [12HT4] , and the lat- 
ter by two surface point defects called boojums [T31 [T7] . 
If the confinement is soft, which happens to be the case if 
the nematic droplet is suspended in a fluid |18j , not only 
the director field adjusts to the presence of the interface, 
but the interface itself adjusts to the director field as well 

Spindle-shaped nematic droplets, or factoids, have in- 
deed been observed in suspensions of elongated colloidal 
particles, under conditions where isotropic and nematic 
phases coexist [19 . Their shape and internal structure 
has been successfully described in terms of the bulk elas- 
tic and surface properties by applying variational theories 
[2] [3l [19] . These make use of plausible families of droplet 
shape and director-field structure, and have been applied 
to extract information on anchoring strengths and elastic 
constants from polarization microscopic images of fac- 
toids ei nana us]- 

Interestingly, values of the ratio of the anchoring 
strength and the surface tension obtained by curve fit- 
ting to experiments were found to be very much larger 



than expected from theoretical prediction [2J El SI] • It 
is not clear if this discrepancy is caused by the varia- 
tional nature of the theories, warranting a more precise 
theoretical investigation of the coupling between direc- 
tor field and droplet shape. In this paper we do just 
that, by invoking Frank elasticity theory and a Rapini- 
Papoular-type description of the interfacial properties of 
the nematic [HI [21]. We restrict ourselves to the two- 
dimensional case, of which physical realizations may be 
found in certain types of Langmuir monolayer [3 1151 - 

UH |23| |2l ESI 123. 

For this kind of two-dimensional system, and in the so- 
called equal constant approximation, we develop a fast, 
and easily implementable numerical method for finding 
the optimal director field for fixed arbitrary shapes and 
arbitrary interfacial energies. Our method takes the form 
of a simple iterative scheme and employs conformal map- 
pings and fast fourier transforms. Hence, the method 
is useful for studies of nematics in confined geometries 
[niEniED] ■ However, we are also able to apply our method 
to study properties of two-dimensional droplet domains 
that are soft and can optimize their shape. Our method 
provides an alternative to numerically more involved fi- 
nite element studies [231 [25] . 

As opposed to earlier quasi-analytical studies [2J El [121 
I20j . which need to variationally optimize both the direc- 
tor fields and domain shapes within certain classes, in 
our case there remains only one parameter to be opti- 
mized, being the shape of the domain. We pick a rich 
class of physically plausible spindle shapes, which con- 
sist of a domain bounded by circle sections, and extend 
it to include smooth, rounded tips. Within this class of 
shape, we compute a state diagram mapping out the di- 
rector field and the domain shape, such as the tip being 
round or sharp, as a function of the various material con- 
stants and domain size. In addition, we show that for the 
spindle shapes with sharp tips, the equilibrium director 
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fields are accurately described by fields constructed from 
circle sections. 

In section [XXJ we first describe the free energy of our 
model and discuss the relevant dimensionless groups de- 



scribing the physics of the problem. Section III discusses 
the equations that determine the optimal directorheld for 
fixed domains, and our light-weight numerical method for 
finding these optimal fields. In section |IV| we put for- 
ward a parametrization of a class of domain shapes that 
interpolate smoothly between sharp and rounded spindle 
domains. In section [V] we present optimal director fields 
calculated for the mentioned class of droplet shapes. Fi- 



nally, in section VI we variationally determine the opti- 
mal shape within our class of domain, and present a state 
diagram of optimal droplet shapes and director fields as 
a function of the relevant material parameters. We end 
the paper with a summary of our findings and a discus- 
sion of possible future extensions of our work in section 



II. FREE ENERGY 

In order to describe the domain shapes and director 
fields, we use a macroscopic theory based on a competi- 
tion between the elastic and surface properties of nematic 
droplets that float freely in an isotropic phase. We model 
this by applying the Frank elasticity theory and a suit- 
able anchoring surface energy that couples the director 
field to the interface, and restrict ourselves to two spatial 
dimensions. 

The director field n indicates the average orientation 
of the particles in the nematic phase and we write it as 



n = (cos 9, sin 6), 



(1) 



where = 0(x, y) is the angle that a particle at position 
[x, y) makes with the horizontal axis, which we choose as 
the axis of mirror symmetry. The free energy F = F e +F s 
of a nematic drop of given surface area S consists of two 
contributions: an clastic contribution F e and a surface 
(boundary) contribution F s : 

i) The Oseen- Frank elastic energy F e is associated 
with the distortion of the director field. In our two- 
dimensional model the deformation energy density com- 
prises only a splay and a bend contribution [21], 



F e = l [ (^i(V-n) 2 + X 3 (V x n) 



dA, 



(2) 



where K\ and K% are the splay and bend elastic con- 
stants, respectively and where the integration is taken 
over the entire surface of the domain O. We adopt the 
equal-constant approximation, where the splay and bend 
elastic constants are presumed to be equal, the sum of 
the elastic deformation energy densities can be written 
as -ftf|VO| 2 /2, with K the average elastic constant. 

it) The interfacial energy F s is associated with the an- 
choring of the director field to the interface (the phase 



boundary). It depends on the angle ipw — between the 
director field 0, and the angle ipN of the surface normal 
with the horizontal axis, and in the most general form it 
is written as 



F. = I a(ip N -e)dS, 
I dn 



(3) 



where <9f2 is the boundary of the domain, and the energy 
density a is given by 

oo 

a(ip N - 0) = a + ^ a " cosn(ip N - 0). (4) 

n=l 

Here, a is the bare surface (line) tension that specifies 
the energy per unit length of the boundary, irrespective 
of the director field 0. The higher order a n specify how 
the boundary energy depends on the director field orien- 
tation. Rudnick and co-workers studied the case of only 
a single a n unequal to zero |15[ I17j and nonzero values 
of both eti and a,2 [24, 25] , If we retain only n = 2 in the 
sum m Eq. @, corresponding to the interfacial energy 
of a two-dimensional nematic [15] . In that case, the in- 
terfacial energy density becomes of the Rapini-Papoular 
type i 



a((fiN - O) = a + a 2 cos 2(tp N — 0), 



(5) 



where a 2 is the anchoring strength of the director field. 
Note that the a2-term can be both positive and negative, 
so in order to prevent unphysical negative total surface 
energy densities, we must have the restriction that \a 2 \ < 
a . 

From now on, we restrict ourselves to the case of planar 
alignment, implying that the director field favours a tan- 
gential orientation to the boundary. However, we shall 
point out that fields with a preference for homeotropic 
anchoring can be easily obtained with our methodology 
as well. 

Combining the above two energy contributions, we 
then have for the total free energy 



F[Q] 



K 



Ve\ 2 dA+ / a(<p N -G)dS. (6) 



no 



We will work with a fixed domain area A$ and all lengths 
are scaled in units of 



P = V 'Aq/tt. 



(7) 



We make the bare surface (line) tension eto and anchoring 
strength a 2 dimensionless by introducing 



and 



7 



2flop 
K ' 



2a 2 p 
K 



(8) 



(9) 



Note that in view of the restriction < a 2 < a we 
always have that 7 < uj. We can write 7 = 2p/£ with 



3 



£ = K/d2 the extrapolation length [2J, so small values 
of 7 should give rise to a homogeneous director field, 
whereas a curved, bipolar field is to be expected for large 
values of 7. 

In order to obtain both the equilibrium shape and di- 
rector field of a domain, one should simultaneously min- 
imize the total free energy Eq. ([6| with respect to O and 
fl. However, this is a formidable mathematical problem 
that is not easily solved analytically [T51 [33], and still 
difficult numerically [24j [25]. Therefore, we first con- 
sider the problem of a fixed shape f2 in the next section, 
and develop an easy-to-implement, numerical method for 
finding the optimal director field that minimizes the free 
energy 



be an analytic function that maps the domain f2 con- 
formally onto the unit disk D : \w\ < 1 and dtt onto 
3D. The process of mapping coordinates z of the do- 
main Q conformally onto coordinates w of the unit disk 
D is depicted schematically in Fig. [I] If, and only if, O 



v> = M 




n d 



III. CALCULATING THE OPTIMAL TEXTURE 
OF A FIXED DOMAIN 



FIG. 1: The conformal map / : C — > C maps a point z £ Q 
bijectively to a point w G D, the unit disk. The inverse / 
performs the opposite mapping. 



If the shape of the domain ft is kept fixed, variational 
minimization of the free energy ([6]) with respect to the 
director field leads to the following boundary value 
problem for O [15J: 



V 2 e(x,y) = 0, 



(x,y) e n 



(10) 



d&(x,y) 
3n 



(x,y)€dn, (11) 



with 



7[©]0> y) = <t'(<Pn(x, y) - Q(x, y)). 



(12) 



Here, <pn is the angle that the normal to the boundary 
makes with the reference axis, and <j'(u) = da/du with a 
the interfacial free energy density (a free energy per unit 
length), Eq. Q. For arbitrary domains Q : the nonlinear 
boundary condition makes the boundary value problem 
posed by Eqs. ( |To] ) and (111 highly intractable for so- 
lution by analytical means. We will therefore resort to 
numerics: using conformal mappings we reduce the prob- 
lem to a single integral equation for the director field O 
on the boundary of the domain, which is solved by an 
iterative procedure employing Fourier transforms. This 
method is fast, easy to implement, and suited for arbi- 
trary domain shapes. In the remainder of this section we 
will outline our solution procedure. 

Inspired by the work of Rudnick and collaborators [H] , 
we start by using the complex variable z — x + iy to rep- 
resent the two dimensional coordinate (x,y) € il, thus 
gaining access to the powerful methods of the field of 
complex analysis. In particular, we use conformal map- 
pings to map the domain f2 onto the unit disk D. A con- 
formal map is an analytic function that maps one domain 
onto another, while preserving harmonicity of functions 
defined on either domain |43j . 

Let 



w = f(z) 



(13) 



is harmonic on S7, then the function 

~M = ©(/-»), weD (14) 

is harmonic on D, where / _1 denotes the inverse opera- 
tion of /. Moreover, if satisfies the boundary condition 
( 11 ), then S satisfies the modified boundary condition on 



3D 



where 



9[ Eo /](/"») ^ m(w)> 



dn D 



(15) 



r)tl 1 



denotes the normal derivative on the unit disk 



D, and the functional g [•](•) is defined in (12). It follows 
that the problem (10 111 on 17 can be mapped onto a 



Neumann problem on the unit disk D: 

( V 2 S = (w eD) 

dE 



. dni 



= h[E](i 



(16) 



(w e 3D) 



with h defined as in Eq. (151. Once the solution 5 to this 



problem is found, the solution to the original problem 
( 10 1, (111 can easily be reconstructed as = 3 o /. 



We solve the above problem by iteration. Starting with 
an initial guess Sq, successive approximations S n _|_i (n = 
0, 1, . . .) to the solution of (16 1 are obtained by solving 
the intermediate Neumann problem 



r v 2 s n+1 = 



(w e D) 



3E 



71+1 



3n D 



h[E n ](w) (we 3D), 



(17) 



where it should be noted that the boundary condition 
is now linear. This problem is an ordinary Neumann 
problem, which is solved by |47) 

E n+1 (w) = -^- f Af(w,w')h[E n ](w , )dw , 1 (18) 

27r JdD 
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where the Neumann function Af(w,w') is specified by 

Af(w,w') =-log(\w-w'\\l-ww'\), (19) 

and where w denotes the complex conjugate of w. 

The nonlinear integral Eq. ( 18 ) provides an itera- 



tion scheme for constructing successive approximations 
of the director field, based only on its values on the 
boundary of the unit disk. It suffices therefore to only cal- 
culate the director field values on the boundary through- 
out the iteration process. Should iteration of the integral 
Eq. ( 18 ) yield a converging sequence 5„ for n — > oo, then 



we have in principle solved problem ( 16 ) 



However, evaluating the integral numerically at each 
iteration step for a discrete set of boundary points turns 
out to be a cumbersome and slow process. Instead, the 
integration step is circumvented by considering a Fourier 
series expansion of h\E n ](w), with w = exp(i^s) € 3D: 



bjj, sin kip. 



(20) 



k=l 



Because of the nematic symmetry of the director field 
and domain f2, which is respected by the conformal map 
/, we only need to consider sinusoidal terms. The Fourier 
coefficients bk, of a function h(ip), can be obtained from 
the forward Fourier sine transform J-fc, which we define 
as 

1 f 2 * 

h = Fk[h{tp)] — - I h(ip) sin kipdip, (21) 
t Jo 

and similarly we define the inverse Fourier sine transform 
J 7 ^ 1 , that reconstructs the function h(ip) from its Fourier 
coefficients bk again: 



Kv) = F<p 1 M = bk sin kip - 



(22) 



fc=i 



Inserting the expansion ( 20 ) in the integral Eq. ( 18 ) , and 
noting that when w, w' £ 3D the Neumann function ( 19 ) 
assumes the particularly simple form 



A/"(e lv , e lv ') = - log (2 - 2cos(v? - </)) . 



(23) 



then the integral in Eq. ( 18 ) can be evaluated analyti- 
cally gO] to yield 



oo , (n) 

' 



Sn+iCO = ]T -j- sinfc^ = F-^/k 



(24) 



k=l 



Using the above calculations and definitions, the suc- 
cessive approximation iteration scheme for the integral 
Eq. (18) turns into an iteration scheme for the expan- 
sion coefficients bk- 



,(«+!) 



7 sin (20 — 2J 7 ^ 1 b^/k' 



\nf- x {^))\- 



(25) 



where is the angle that the normal to 3ft makes with 
respect to the x-axis, evaluated at the point f~ 1 (e 1(p ). 
Finding the optimal director field has now been reduced 
to an iteration of the above equation, and when the se- 
quence converges the director field can simply be found 
by first calculating H through Eq. (24), and subse- 
quently evaluating = So/. Each iteration step in 
the above scheme requires only two Fourier transforms, 
which can be done both accurately and fast using Fast 
Fourier Transforms (FFTs) [41]. Moreover, the above it- 
eration scheme is easy to implement, making the determi- 
nation of the optimal director field for a fixed shape with 
a known conformal map a simple computational task. 

We conclude this section with some notes on practi- 
cal issues. Firstly, the iteration scheme (25) is adjusted 
slightly in order to aid convergence. 



ing an approximation b 



hand side of Eq. (25) 



Rather than replac- 
in its entirety by the right- 



the new approximation u k 



(n+1) 



is then composed as a combination of /i times the right 
hand side of Eq. (25 1 and (/i — 1) times the previous ap- 

The 



proximation b k n \ with some parameter < fx < 1. 
convergence of this sequence is then accelerated using a 
Shanks transformation [42] . which roughly doubles the 
speed of convergence. 

Secondly, when we obtained a solution S„ + i to the in- 



termediate Neumann problem (17) at iteration step n, we 
have in fact obtained a whole family of solutions. Addi- 
tional solutions can be found by simply adding a constant 
c G [0, 27r], for S„+i -I- c still solves problem (17). The 



constant c determines the average value of the director 
field taken over the entire domain. However, when the 
sequence S„ has converged to some final value S^, the 



same does not hold true for the original problem (16). 
For if Soo solves (16), + c does in general not solve 



that same problem due to the nonlinearity of the bound- 
ary condition. But problem ( 16 ) in fact does allow for 
multiple solutions. Solution ( 18 1 represents a director 



field which has an average orientation of (S) = 0, so, in 
other words, we have tacitly picked the c = solutions. 

To find director fields S with a different average orien- 
tation (3) = c, we need to select the corresponding solu- 
tion with offset c at each and every step of the iteration. 
For domains f2 with symmetry along the x and y axis, the 
only two possible values of c are c = and c = 7r/2 (mod- 
ulo it/2). This can be seen from the Neumann condition 
f an h[E](w)dw — that has to be satisfied for problem 
( 16 1 to have a solution [44]. For elongated domains along 
the x-axis it turns out that with a-i > the c = solu- 
tion corresponds to an energy minimum and c = 7r/2 to 
a maximum. When a-2 < the opposite is true and we 
should select the c = ir/2 solution at each iteration step. 

This concludes the description of our methodology. 
Note that this method is in principle applicable to any 
domain f2, and any arbitrary boundary energy functional 
g[Q] (z). However, in this paper we shall restrict ourselves 
to nematic phases, with boundary energies of the form of 
Eq. Q. To proceed further, we now need to specify the 
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FIG. 2: Spindle shaped domain with with aspect ratio e s = 2. 
Relevant variables in the parametrisation and conformal map 
are indicated. The symbols zi, .., Za that are drawn at a point 
(a;, y), correspond to complex numbers z = x + iy. 



exact shape of the domain Q. 



IV. DOMAIN SHAPES 

We introduce two classes of domain, as well as their 
parametrization and conformal mapping, for which we 
will calculate the optimal director field as outlined in Sec. 
Ill In the next sections, we use these classes of domain 
to approximate the optimal shape. The first class of do- 
main is that of spindle-shaped domains with sharp tips, 
whereas the second class contains domains with rounded 
tips, and these are constructed from the first class, to be 
discussed below. 

The choice for the spindle shapes is inspired by ex- 
perimentally observed domain shapes in two 
and three [45J 06] dimensions, and the optimal shapes 
found in numerical studies [53J [25J . Sharp-tipped spin- 
dle shapes are domains bounded by two identical circular 
arcs, as depicted in Fig. [2j The circular arcs making up 
the spindle shape have a radius Rq, and intersect each 
other at the point z — ±Z\ (recalling that z = x + iy), 
located on the x-axis. The y-axis is intersected at point 
±z 2 , and provided that \z%\ < \z\\, we can define the as- 
pect ratio of the spindle domain as e s = |zi|/|z2|- The 
spindle shape in Fig. [2] has an aspect ratio of e s = 2. 
Using the area A$ of the domain and aspect ratio e s as 
input parameters, the spindle boundary is parameterized 
in polar coordinates, in the first quadrant, by 

R(ip) = Ro(e s ) (V 1 - 3 2 cos 2 <p - 8 sin (26) 

where tp is the polar angle, and 

e 2 - 1 

3 = 

' e 2 + 1 

The circular arc radius is given by 



(27) 



Ro(e s ) = 



A. 



o 



^ 2 (arccos 3 - B^J\ - /3 2 ) 



(28) 



The normal angle to the surface at a point (R(ip),tp) is 
equal to 



ifN = p — arctan 



1 OR 



(29) 



The above parametrizations are restricted to the first 
quadrant < ip < ir/2, but other quadrants can be found 
by simply mirroring the first quadrant. 

The conformal mapping of the spindle shaped do- 
main can be constructed from two consecutive mappings. 
First, the spindle domain f2 is mapped to the upper half 
plane u 6 C, Im(u) > by the conformal map [S3] 



(z - Zx)(za - z 3 ) 



7T /2a 



zen, (30) 

R(Q), z 2 = iR(tt/2), 



u ( z ) = -i/ 

V (z ~ Z 3 )(Zl - z 3 ) 

with the complex numbers z\ 

z 3 = — z\ (see Fig. [2]), and the locus of the upper circle 
section 2:4 = iR(tt/2) — iRo(e s ). The angle 2a is that of 
the sharp tip of the spindle, as indicated in Fig. [2] and 
is calculated as 

a = arcsin ( — ) = arccos 3, (31) 
\Ro(£s)J 



with 8 as defined in Eq. (27). Next, defining 



and 



v(u) = 



9{u) 



U — U\ W3 — Wl 

u - u 2 w 3 — w 2 ' 



Wi — v(u)iV2 
1 — v(u) 



(32) 



(33) 



then w = g (it) is the conformal map that maps the upper 
half plane to the unit disk in the w-plane, where the 
points u\ , u 2 and 00 of a line in the upper half plane are 
mapped onto the points w\ = 1, 102 = 1 and w 3 = —1 on 
the unit disk. When we pick u% = u(z±) and u 2 = u{z 2 ), 
then the composite map (see also Fig. [I]) 



f{z)=g{u{z)) 



(34) 



maps the spindle domain onto the unit disk, where the i- 
th quadrant in the z-plane is mapped to the i-th quadrant 
in the to-plane (i = 1, 2, 3, 4). 

The inverse map, mapping the unit disk conformally 
to a spindle domain, can be computed by inverting Eq. 
(34), resulting in 



(\ 2a I tt 



1 



(35) 



-1 



Interestingly, through the introduction of an additional 
parameter A, with < A < 1, a whole new family of 
domain shapes Q\ can be generated with the conformal 
mapping 



zx = A x f-\\w), 



(36) 
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FIG. 3: The family of s mooth spindle shapes generated by 
the conformal map ( |36| l, for various values of A between 
and 1 with outward lying shapes for increasing A (indicated 
by arrow). The outermost shape drawn with a dashed line, 
is the original, sharp-tipped spindle shape with aspect ratio 
e s = 2 and A = 1. 



which maps only a part of D onto the z-domain, creating 
a more rounded spindle shape without a sharp tip. The 
new smoothed spindle shapes are generally smaller than 
the original spindle shape, and a constant A\ is added 
for area conservation. This class of shapes transforms 
fluently from the sharp-tipped spindle domain (A = 1) 
to the unit disk (A — > 0), and they bear strong resem- 
blance to the fully numerically calculated optimal results 
[24"t 125] . The family of smooth spindle shapes is shown 
in Figure [3j where we have temporarily set A\ = 1 for 
clarity. 

The proper value of A\ can be calculated by noting 
that the droplet area A can be calculated using the di- 
vergence theorem: 

2A = f V-(x,y)dA= [ (x,y)-ndS = Alh, (37) 
Jn x Jdn x 

with 

h = J (Rc ^(0)lm ^ - Im V(#)Re ^ ] d0, (38) 



89 I 



and where we have written tp(0) = / _1 (Aexp(i0)). In 
order to conserve the area of the droplet Aq when varying 
A, we thus must have that 



and the radius of curvature R c at the tip of the domain, 

%KA\\e r a 



R c {£si A) 



tt(1 + K) 2 (l + A) 2 + 4aX(K 2 - 1) ' 



(41) 



with 2a the tip angle as in Eq. ( |3l[ ), and where we have 
defined 



K = 



e, — l \ /A-l x 



e s + ij V A + 1 



(42) 



The radius of curvature as defined above is defined in 
units of the waist size |/ _1 (iA7r/2)| of the droplet, such 
that it is always a number between and 1. Also, in the 
case of A = 1, we have that e r = e s , and in the remainder 
of this paper we shall use e r to indicate the aspect ratio 
of both round-tipped and sharp-tipped shapes. 

Having introduced the class of shapes of interest, it be- 
comes important to note an implementation issue for the 
solution method of Sec. |III| for elongated domains ft with 
sharp tips, such as spindle shapes with high aspect ra- 
tios. Conformal mappings between an elongated domain 
and the unit disk suffer from a phenomenon called crowd- 
ing [SD1 [HI] • The concept of crowding can be understood 
as follows. Consider a set of uniformly distributed points 
zgSl. When fl is highly elongated, the density of images 
w = f(z) of these points will become highly non-uniform 
in the unit disk. Conversely, a uniformly spaced grid 
on the unit disk gets mapped to a highly non-uniform 
spaced grid in the domain J7. The number of gridpoints 
on the unit disk required to maintain a certain density of 
images on scales exponentially with the aspect ratio of 
the domain [50] . 

Additional crowding originates from highly distorted 
boundary curves dfl |52j . as occurs for instance in the 
case of a sharp tip in the domain. For the spindle 
shaped domains discussed in the next section, our so- 
lution method, which requires equidistant grid points on 
dD on behalf of the FFT, needs in the order of up to 10 6 
gridpoints to correctly sample the director field in the 
tip of the domain when the aspect ratio of the droplet 
domain is valued around 3. Due to the efficiency of the 
FFTs such large grids do not pose a problem in terms 
of computation time, but the exponential scaling of the 
number of required gridpoints does prevent us from ac- 
curately investigating higher aspect ratios than 3. 



Ay 



2Ao 
h ' 



(39) 



The parameters e s and A provide not much insight 
into the physical shape of the domain. Moreover, the 
droplet shape is very (in)sensitive to small changes in A 
for A close to 1(0). Therefore, we introduce two new and 
more convenient descriptors, being the aspect ratio of the 
round-tipped droplet 



e r — e r (e s , A) 



r x (A) 



|/-HiA7r/2)| : 



(40) 



V. OPTIMAL DIRECTOR FIELD FOR FIXED 
SPINDLE SHAPES 

In this section we employ our solution procedure of 
Sec. |III[ to calculate the optimal director field for a fixed 
shape from the spindle classes introduced in the previous 
section. Before turning to the general case, it is inter- 
esting to consider the simplest spindle shape: a circu- 
lar domain. For this domain the conformal mapping is 
simply the identity operation. Despite the complicated 



nonlinear boundary condition in Eq. (12 1, Rudnick and 
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Bruinsma |15) found an exact analytical solution on the 
unit disk D, provided do and only one other in the 
expansion Q is nonzero. Their solution is of the form 



©0, y) = -(q(x + w) - q(x - w)), 



with 



q(z) = - log 1 



and the length L given by 



z 



L 2 = 1 + 

7 



(43) 



(44) 



(45) 



In the case of a nematic liquid crystalline phase, for 
which k = 2, it can be shown that the above director field 
is equal to a bicircular field that has field lines consisting 
of circular sections intersecting the x-axis at points ±L, 
as illustrated in Fig. |4j These virtual defects lie outside 
the domain Q, and are commonly referred to as boojums 
in the limiting case for which they reside exactly on the 
boundary of the domain. 

The representation of such a bicircular field in polar 
coordinates (r, tp) is given by 



0(r, tp) = arctan 



-2r 2 



cos tp sin tp 



r2 (sin tp 



cos 2 tp) 



(46) 



As we shall see, this field is a very good approximation to 



the solution of the problem ( 10 ), ( 11 ) if the domain f2 is 



and a field that is parallel to the boundary. For more 
extreme values of the parameter 7, we find that for 7 — > 
the director field becomes homogeneous, whereas for 
7^1 the director field becomes rigidly anchored to the 
boundary. 

This behaviour is best observed by evaluating the di- 
rector field at the boundary of the domain, as shown in 
Fig. [6j Here, the value of the director field on the bound- 
ary of the domain Q is plotted as a function of the polar 
angle ip in the first quadrant, and for various values of 7 
ranging between 7 = 0.1 and 7 = 10. The (red) dashed 



spindle shaped, to which we shall now turn our attention 





FIG. 4: The field lines of a bicircular field are circle sections, 
intersecting the i-axis at points ±L. Bicircular fields corre- 
spond to the lowest energy configuration of the director field 
for circular droplet domains. 

Starting with spindle shapes with sharp tips, a typical 
example solution of the problem (10 1, ( flT| ), is shown in 
Fig. [5]j a) by means of a plot of director field lines. Here, 
the anisotropic line tension (anchoring strength) is set 
equal to 7 = 2 [see Eq. ( 9 1] , and the domain aspect ratio 
fixed at e s = 2. Figure 5 b) shows the optimal director 
field for the same parameters, but with a round-tipped 
shape. For these parameters, the director field seeks a 
compromise between an undeformed, homogeneous field, 



FIG. 5: Field lines of the optimal director field for 7 = 2 
and a fixed shape within the class of smooth spindle shapes: 
(a) sharp-tipped shape with aspect ratio e r = 2 and R c — 0, 
and (b) round-tipped shape with e r = 1.85 and R c = 0.21, 
corresponding to e s = 2, and A = 0.99. 

line in Fig. [6] indicates the limiting case of a field parallel 
to the boundary at all points. For sharp-tipped shapes, 
at tp — 0, this field goes to the value a as defined in Eq. 
(31 1 (see Fig. [6]ja)), while for round tips the parallel field 
always goes to —it/2. 

The optimal fields for sharp-tipped shapes turn out 
to be very similar to the bicircular director field from 
Fig. [4] Indeed, when we restrict the director field to 
be bicircular and minimize the free energy with respect 
to position L of the (virtual) boojum, we find that the 
optimal bicircular field is numerically indistinguishable 
from the calculated optimal director field. The relative 
differences in free energy between the optimal field and 
the bicircular field are shown in Figure ^ a), for a large 
range of values of 7 and aspect ratios s s . There is a 
maximum at 7 ~ 0.5 and e s ~ 1.3, where the relative 
energy difference is ~ 3.5 • 10 -4 . For all other parameter 
values the bicircular field performs even better. 

However, the boundary condition specified by Eq. ( 11 ) 




FIG. 6: The angle <d of the optimal director field with the 
long axis of the droplet, evaluated at the boundary of the do- 
main, and plotted as a function of the polar angle for various 
values of 7 ranging between 0.1 and 10. The arrow indicates 
the direction of increasing 7. Dashed (red) line indicates a 
director field perfectly aligned to the boundary, as would oc- 
cur in the limit 7 — > 00. Shape parameters are identical to 
those of Fig. [5] with (a) a sharp-tipped shape with aspect 
ratio e r — 2 and R c — 0, and (b) a round-tipped shape with 
e r = 1.85 and R c = 0.21. 
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FIG. 7: Relative energy difference AF = (F c - F)/F, of the 
optimal director field with energy F and the optimal bicircu- 
lar field with energy F c , shown as a function of the domain 
aspect ratio e a , e r , and the line tension anisotropy 7. Panel 
(a): sharp-tipped spindle shapes, (b): round-tipped spindle 
shapes with a fixed radius of curvature R c = 0.2. 



cannot be satisfied for all values of <p for a single choice 
of L. This can be seen by solving for L, for particular 
choices of <p for which the equations simplify, e.g., ip — > 
and if = 7r/4. For these two cases, two different results 
for L are obtained. Thus, although the likeness of the 
two fields is tantalizing, this proves that the bicircular 
field is not an exact analytical solution to the boundary 
value problem (10 11 1, on a spindle shaped domain. For 



all practical purposes however, the usage of the spindle 
field instead of the optimal solution appears justified for 
spindle shapes with sharp tips. The reduction of the field 



optimization to a single parameter to be optimized for 
makes this especially convenient. 

For domains with round tips the bicircular field is not 
nearly such a good approximation to the optimal field 
as it is for sharp-tipped domains. This is evidenced in 
Fig. [7][b), where the relative energy difference between 
the optimal field and the bicircular field is plotted against 
7 and the aspect ratio e r , while maintaining a constant 
radius of curvature R c = 0.2. For large values of 7, the 
bicircular field cannot follow the domain boundary in the 
tip very well, leading to relatively high energy differences. 
But even for the rounded tips, there exist large regions in 
parameter space where the bicircular field approximates 
the optimal field very well. For smaller radii of curvature 
R c , the bicircular field performs better, whereas for larger 
R c it performs worse. 

Having established that the bicircular field approxi- 
mates the optimal field very well for sharp-tipped shapes, 
we can characterize the optimal field by means of the vir- 
tual boojum locations L of the bicircular field. Figure [8] 
shows a plot of the optimal value of the location L of 
the (virtual) boojum, as a function of the shape eccen- 
tricity e s and various values of the anisotropic boundary 
energy strength 7. At e s = 1 L obeys Eq. (45). Clearly, 



the higher 7, the closer the virtual boojums are located 
to the tip of the droplet, until in the limit 7 — > 00 the 
boojums are located exactly on the tips of the spindle 
domain. 
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FIG. 8: Virtual boojum location L of the optima l bi circular 
field, in units of the location R(<p = 0) [see Eq. (26l] of the 
tip of the droplet, and as a function of droplet aspect ratio 
e s , for logarithmically spaced values of 7 between 0.1 and 10. 
The arrow indicates the direction of increasing 7, with the 
7 = 1 case indicated by a dashed line. 



VI. THE OPTIMAL SHAPE 

In this section we determine the optimal shape within 
the class of sharp and smooth spindle shapes by varia- 
tionally minimizing the free energy, Eq. ([6| , with respect 
to the shape parameters e r and R c given by Eqs. (40) 
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and (41 ). At each step in the minimization procedure, we 



III to constrain the 



employ our solution method of Sec. 
director field to the optimal field for the current shape 
parameters. 

Results of this optimization procedure are given in Fig- 
ure |9j showing the optimal shape for various values of uj 
between 0.01 and 100, and < 7 < uj. Each line corre- 
sponds to a particular value of uj, where arrows are used 
to indicate the direction of increasing uj, and increasing 7 
(for Fig. [9](d)) . Panel (a) shows the radius of curvature 
of the tip of the domain, as a function of the ratio j/uj, 
which measures the anisotropy of the interfacial energy. 
For large enough 7/w, R c becomes nought, indicating a 
transition to a sharp tip. Panel (b) shows the tip angle 
of the optimal shape for the case of sharp-tipped shapes 
indicated with dashed lines. The solid lines in this panel 
are extrapolations, and indicate tip angles if the tips were 
forced to be sharp, i.e., if we ignore the possibility of the 
rounded tip. Panel 9]jc) shows the aspect ratio of the op- 
timal shape, with the convention that solid lines indicate 
shapes with rounded tips, using the same data points as 
those in panel (a), whereas dashed lines indicate sharp- 
tipped shapes corresponding to the data points shown in 
panel (b). 

These figures show that for 7 — » 0, the boundary en- 
ergy is completely determined by the isotropic line ten- 
sion uj and the domain shape becomes perfectly circular 
in order to minimize its circumference, as evidenced by 
the radius of curvature R c — > 1 and e r — > 1. As 7 is 
increased, the tip becomes sharper (decreasing R c ) while 
simultaneously the shape elongates (increasing e r ), until 
at some critical value of 7 < uj the radius of curvature 
R c = 0, and the shape tip becomes sharp. It should be 
noted that the transition point between round and sharp 
tips occurs at different values of j/uj when uj is varied. 
In the limit uj —¥ 00, the transition occurs immediately 
at 7/w = (the optimal shape is always sharp-tipped), 
and the transition point moves up monotonically towards 
7/0; ~ 0.45. Continuing with increasing 7 beyond the 
transition point, the tip remains sharp and the shape 
will continue elongating, until in the limit 7 — > uj the 
aspect ratio e r — > 00. 

Finally, the two shape parameters for round-tipped 
shapes are plotted against each other in the inset Fig. 
[9jd) , showing clearly the tendency of increasing 7 to elon- 
gate the droplet and sharpen the tip, whereas increasing 
uj tends to favour more circular domain shapes. Addi- 
tionally, the inset contains no data points in the upper 
right part of the figure. The corresponding shapes would 
have high aspect ratios but a radius of curvature compa- 
rable to the droplet waist (R c ~ 1), and apparently there 
exist no values of 7 and uj for which such shapes become 
energetically favorable. 

The results of our numerical calculations of the opti- 
mal director-field structure and shape (within our class 
of smooth spindle shapes), can be gathered in a single 
phase diagram as shown in Figure [TO} We distinguish six 
regimes, labeled I through VI, based on three defining 




y / co 

FIG. 9: The optimal domain shape, characterized by its as- 
pect ratio, e r , and radius of curvature at the tip, R c , plotted 
for various values of the dimensionless isotropic line tension 
uj — 0.01 . . . 100, logarithmically spaced with the direction of 
increasing uu indicated with an arrow. Panel (a) shows R c as 
a function of the ratio of anisotropic and isotropic line ten- 
sion 7/w, and panel (b) the tip angle of the optimal shape, as 
a function of 7/cj, indicated with dashed lines. Solid lines in 
this panel are extrapolations that indicate the tip angle where 
the shape constrained to have a sharp tip. Panel (c) shows 
the aspect ratio s r of the optimal shape as a function of 7/w, 
where solid and dashed lines are used for round and sharp 
tips, respectively. In panel (d) R c and e r are plotted against 
each other, the arrows indicating the direction of increasing 
material constants oj and 7. 



properties: 

i) The director field being homogeneous or curved. For 
small values of 7, the bulk elasticity dominates over the 
interfacial anchoring, so the director field is rigid and ho- 
mogeneous. For large values of 7, the anisotropic bound- 
ary energy dominates and the director field will accom- 
modate by aligning itself parallel to the shape boundary. 
The (smooth) transition between the two types of field 
is chosen to occur at 7 = 1, and is indicated by the solid 
(red) line in Fig. 
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Points to the left of this line cor- 
respond to homogeneous director fields, whereas points 
to the right have a bipolar (curved) director field. The 
crossover value 7 = 1 is sensible because the anchoring 
energy and the elastic deformation energy are then ap- 
proximately equal. 

ii) The droplet tip being rounded or sharp. For small 
values of the dimensionless anisotropic line tension 7, the 
dimensionless isotropic line tension uj dominates the total 
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FIG. 10: The shape and director field of the nematic do- 
mains are shown as a function of the dimensionless anchor- 
ing strength j/uj and dimensionless surface tension uj. Six 
regions, indicated by the Roman numerals, are found that 
follow from three characteristics we distinguish. The dotted 
line gives the boundary between a smooth and a sharp tip 
of the domain; the dashed line shows the transition from an 
elongated to a round domain; the solid line is the boundary 
between a uniform and a bipolar director field. 



interfacial (line) free energy. The shorter total boundary 
length associated with rounded tips outweighs the rise in 
anisotropic surface energy due to the larger attack angles 
at such a round tip. Hence, for small values of 7 the 
droplet tip becomes rounded, whereas for higher values 
of 7 the opposite situation occurs and the tip becomes 
sharp. In Fig. [10] the transition from round to sharp is 
indicated by a black, dotted line. Points above this line 
have sharp tips, points below have rounded tips. 

Hi) The elongation of the droplet. The elongation of 
a droplet is defined by its aspect ratio e r . [Note that 
e r = s s for sharp-tipped shapes.] For high isotropic line 
tensions u, the droplets will tend to be close to circular to 
minimize their interfacial energy. If the anisotropic line 
tension 7 increases and approaches w, the director field 
aligns along the boundary. If the isotropic line tension 
is small, the director field is rigid. Aligning the field by 
elongating the droplet becomes energetically favourable 
over deforming the director field, at the expense of in- 
creasing the boundary length and associated isotropic in- 
terfacial energy. In Fig. 10 the dashed (blue) line marks 



points with constant aspect ratio, in this case equal to 
e r = 2. All points located above this line have larger 
aspect ratios and are more elongated; all points below it 
have lower aspect ratios and are more round. 

In principle, these three properties allow for eight dis- 
tinct topological combinations in the phase diagram. 
However, we only distinguish six possible regimes in the 
phase diagram, defined by the limiting behavior of the 
dividing lines. These six regimes extend all the way to 
the extreme values of j/lj, and logo; — > ±00. These 



six regimes are labeled I through VI in fig. [lOj and are 
defined by the following characteristics: I) elongated do- 
main, with sharp tip and homogeneous director field, II) 
elongated domain, with sharp tip and curved director 
field, III) round domain, with sharp tip, and curved direc- 
tor field, IV) round domain, with round tip, and curved 
director field, V) round domain, with round tip, and ho- 
mogeneous director field, and finally, VI) elongated do- 
main, with round tip, and homogeneous director field. 

The remaining two topological combinations are elon- 
gated round-tipped domains with a curved director field, 
and round droplets with a sharp tip and homogeneous 
director field. These regimes would constitute a finite 
domain in the center of the phase diagram. However, the 
lines defining transitions between homogeneous or curved 
director fields (i), and round and elongated shapes (iii), 
are somewhat arbitrary as these transitions are smooth 
and the location of the transition is a matter of defini- 
tion only. By choosing different values for the critical 
parameters, these lines can be moved and the phase dia- 
gram will look differently in the center region. Note that 
the transition between round and sharp tipped shapes is 
sharp and resembles a second order phase transition. 

The current critical values are chosen such that the 
dividing lines all cross in a single point. Moving one of 
the lines will immediately introduce one of the remaining 
two topological combinations. But these regions typically 
would only occupy a small area of the phase diagram and 
are subject to arbitrary criteria. By picking the current 
values, the phase diagram only shows the essential topo- 
logical combinations that would be always present as a 
limiting case for any reasonable choice of critical values. 



If we compare the state diagram Fig. 10 for two dimen- 
sional nematic domains with that calculated by Prinsen 
et al. [3] for three dimensional nematic drops, we find 
that they are remarkably similar. Our phase diagram 
has one extra attribute not present in that of Prinsen et 
al. [5], being the transition from round to sharp tipped 
shapes. This generates two more regimes not present in 
the earlier work: regimes IV and VI, i.e., rounded shapes 
with curved director fields, and elongated domains with 
round tips and uniform director field. Interestingly, the 
crossover line for the transition between elongated and 
round shapes coincides to within numerical precision, 
even within regime VI. 

An important point to note is that for a homogeneous 
director field, i.e., in the limit uj — > 0, the optimal shape 
can be obtained from the so-called Wulff construction 
|33j . This shape is then completely determined by the 
anisotropy j/uj of the interfacial tension. However, the 
Wulff shape is not part of our class of spindle shapes and, 
in particular, the limiting value of the sharp-round tran- 
sition in Fig. 10 of j/uj ~ 0.45 does not occur at the value 



of j/uj = 1/3 as predicted by the Wulff construction. 

The fact that the Wulff shape is not included in the 
spindle shape class can be seen as follows: the Wulff 
shape found just above the transition point has an as- 
pect ratio of 2, while the tip angle a wu iff — > ir/2. A 
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sharp-tipped spindle shape of aspect ratio 2 would neces- 
sarily have a tip angle a ~ 0.3-7T, as dictated by Eq. (31 1 



However, the spindle shapes we find for a uniform direc- 
tor field have a free energy that is at most 0.1% larger 
than that of the Wulff shape (at least for the equivalent 
shapes in three dimensions [2]). In the discussion in the 
next section we will elaborate on how the smooth spindle 
class can be expanded to include shapes with arbitrary 
aspect ratios and tip angles, and approximate the Wulff 
shape arbitrarily closely. 



VII. CONCLUSIONS AND DISCUSSION 

In this paper we have presented a general approach 
to find the optimal shape and director field of two- 
dimensional nematic domains. Using conformal map- 
pings, we can transform the original boundary value 
problem with nonlinear boundary conditions on the gen- 
eral domain f2, to a boundary value problem on the unit 
disk. Next, by noting that the value of the director field 
on the boundary uniquely determines the solution, we 
can rewrite the problem in terms of the Fourier coeffi- 
cients of the value of the director field on the boundary 
of the unit disk. Finally, we numerically solve for these 
coefficients with an iterative procedure, which is easily 
implemented and can be done at relatively low computa- 
tional cost due to the usage of Fast Fourier Transforms. 

Using this methodology, we investigate the optimal 
field within a realistic class of spindle shapes, which can 
have sharp as well as round tips. We find that the optimal 
director field for spindle-shaped domains with a sharp tip 
is numerically equivalent to a field that consists of circle 
sections. This means that for practical purposes such a 
bicircular field can be used to model the director field 
of spindle shaped domains of particles that prefer planar 
alignment of the director field to the boundary. 

By establishing the connection of the optimal direc- 
tor field to a given shape, we can search for the optimal 
shape within our class of smooth spindle shapes for given 
material constants. The results for the optimal shape are 



compiled into a single phase diagram (Fig. 10), in which 
we distinguish six characteristic regimes of domain shape 
and director field. The optimal shape and fields are clas- 
sified based on whether the domains are elongated or 
round, have a sharp or a smooth tip, and have a bipo- 
lar or a homogeneous director field. We find that if the 
anisotropy of the interfacial energy is small, the domain 
has a round tip. 

Although we have limited ourselves in this work to 
realistic shapes based on domains freely floating in an 
isotropic background, our approach is straightforwardly 
extended to cover several other interesting related sys- 
tems. For instance, our methodology is applicable to 
any shape for which the conformal mapping to a circle is 
known. A possible application of this is, e.g., the study 
of nematics in confined geometries, or the coalescence 
of two drops by taking a peanut-shaped domain. This 



latter domain can for instance easily be constructed by 
a conformal map of the form f(z) cx tan(A.z), with A a 
parameter between and 1. Alternatively, the class of 
smooth spindle shapes can be easily extended to a more 
general class. 

We have only considered a preferred planar alignment 
of the director field to the domain surface. However, the 
preferred type of anchoring could be perpendicular to the 
boundary (i.e., homeotropic) . In that case 02 < in our 
boundary energy density Eq. ^ , and typically we expect 
a radial director field with a hedgehog defect at the center 
of the domain, at least for large enough domains [TJ] . 

The defect free case with homeotropic anchoring, i.e., 
for small enough domains, is readily covered by our 
method as briefly explained at the end of section |III| 
More generally speaking, our calculations make no ex- 
plicit use of the functional form of the boundary energy. 
They could therefore identically be applied to any other 
type of boundary-energy prescription. 

To account for defects in the interior of the domain, our 
method can be extended because there exist conformal 
maps that map between the unit disk and a disk with a 
small hole cut out at an arbitrary location [S3]. These 
could then be used to generate domains with vanishingly 
small holes, and by setting Dirichlet boundary conditions 
at the boundary of such a hole it would be possible to 
emulate the presence of a defect in the field. 

Unfortunately, there are a number of caveats associ- 
ated with our methodology. Firstly, the crowding phe- 
nomenon of the conformal mapping, as discussed at the 
end of Sec. |IV[ makes it computationally difficult to ac- 
curately calculate the director field for domains with high 
aspect ratios larger than approximately 3, and with very 
sharp tips. 

Secondly, we found that the optimal shape for homoge- 
neous director fields, as given by the Wulff construction, 
possesses particular combinations of aspect ratio and tip 
angle that are not supported by our smooth spindle shape 
class. However, this caveat can be circumvented by se- 
quentially applying more than one conformal map of a 
unit disk to a spindle domain. This way, it is possible to 
generate spindle-like domains, but with arbitrary aspect 
ratios and arbitrary tip angles, including smooth tips. 

Thirdly and finally, the method as presented in this pa- 
per works only in the equal-constant approximation, be- 
cause only then the boundary value problem defining the 
director field simplifies to a Neumann problem. Likewise, 
similar issues arise in generalisations to three dimensions, 
even if we invoke cylindrical symmetry. This means that 
for these cases the method cannot be directly applied, 
but needs to be modified. We leave this for future work. 
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